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ABSTRACT 

We present the first detailed and homogeneous analysis of the s-element con- 
tent in Galactic carbon stars of N-type. Abundances of Sr,Y, Zr (low-mass s- 
elements, or Is) and of Ba, La, Nd, Sm and Ce (high-mass s-elements, hs) are 
derived using the spectral synthesis technique from high-resolution spectra. The 
N-stars analyzed are of nearly solar metallicity and show moderate s-element en- 
hancements, similar to those found in S stars, but smaller than those found in the 
only previous similar study (Utsumi 1985), and also smaller than those found in 
supergiant post-AGB stars. This is in agreement with the present understanding 
of the envelope s-element enrichment in giant stars, which is increasing along 
the spectral sequence M— ^MS^S^SC^C during the AGE phase. We compare 
the observational data with recent s-process nucleosynthesis models for different 
metaUicities and stellar masses. Good agreement is obtained between low mass 
AGB star models (M < SM©) and s-elements observations. In low mass AGB 
stars, the ^^C(q;, nY^O reaction is the main source of neutrons for the s-process; 
a moderate spread, however, must exist in the abundance of ^^C that is burnt 
in different stars. By combining information deriving from the detection of Tc, 
the infrared colours and the theoretical relations between stellar mass, metallic- 
ity and the final C/0 ratio, we conclude that most (or maybe all) of the N-stars 
studied in this work are intrinsic, thermally-pulsing AGB stars; their abundances 
are the consequence of the operation of third dredge-up and are not to be ascribed 
to mass transfer in binary systems. 

Subject headings: nucleosynthesis — stars: abundances — stars: AGB — stars: 
carbon 



1. Introduction 

It is known that the chemical composition of the interstellar medium is oxygen-rich. 
Hence, the overwhelming majority of the stars are formed with a carbon to oxygen ratio 
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lower than unity, and most of them do not change this property during their evolution. 
However, there are exceptions: several classes of stars are known, whose carbon to oxygen 
ratio in the envelope is larger than unity (by number of atoms). These stars are named 
carbon (C) stars. The existence of carbon stars must be related to some specific mechanism, 
acting on a limited population of objects. There are basically three possibilities: i) the car- 
bon enrichment is the result of a deep mixing process suitable to pollute the photosphere 
with the carbon synthesized by shell He-burning in the stellar interior; ii) it is due to mass 
transfer of carbon-rich material after the stellar birth; iii) it dates back to the star's birth, 
due to an anomalous composition of the parental cloud, more enriched in carbon than in 
oxygen. Concerning the last hypothesis, so far no such interstellar clouds have been detected 
in Galactic disk environments, though the hypothesis is not completely ruled out for certain 
low mctallicity stars (Bcvcridgc & Sncden 1994). In the second case, the origin of the carbon- 
rich material is simply moved to another place, i.e. the carbon-rich primary component of a 
binary system. For instance, the stars of the CH spectral type (Luck & Bond 1991; Vanture 
1992a,b,c) most likely owe their carbon enhancement to the transfer of carbon-rich material 
from a companion (now a white dwarf). These stars are usually named extrinsic asymptotic 
giant branch (AGB) stars to distinguish them from those deriving their carbon enhancement 
from nuclear processing in their interiors (hypothesis i), which are called intrinsic. Here, 
we shall deal with a specific subclass of carbon stars, those of spectral type N. We shall 
see later that most probably all of them can be classified as intrinsic C-rich objects. They 
are formed through the mixing into the envelope of newly produced ^^C from He burning. 
Stellar evolution then limits the candidates to evolved stars with masses 1 < M/Mq < 8, 
in particular those ascending AGB. The structure of an AGB giant is characterized by a 
degenerate CO core, by two shells (of H and He) burning alternatively, and by an extended 
convective envelope. In the HR diagram, these stars lay close to the brightest part of the 
Hayashi line. They become long period variable of the irregular, semiregular or Mira types, 
presenting large mass-loss rates: 10~^ to 10~^ Mg/yr (Wallerstein & Knapp 1998). As a con- 
sequence, a thick circumstellar envelope eventually forms, sometimes developing detached 
shells. Depending on its chemical composition and optical thickness, this circumstellar ma- 
terial can obscure partially or completely the central star at optical wavelengths (Knapp & 
Morris 1985; Olofsson et al. 1993; Marengo et al. 2001). Schwarzschild & Harm (1965) 
early showed that thermal instabilities in the He-shell (thermal pulses, TP) occur periodi- 
cally during the advanced phases of AGB evolution. During a TP the whole region between 
the H shell and the He shell (called "the He intershell") becomes convective. After each 
TP, the convective envelope penetrates downward dredging-up material previously exposed 
to incomplete He-burning conditions. This phenomenon is called third dredge up (TDU), 
and its main consequence is the increase of the carbon content in the envelope so that, 
eventually, the C/0 ratio can exceed unity and the star becomes a carbon star. In such a 
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way, the carbon content in the envelope is expected to increase along the spectral sequence 
M^MS— >S— >^SC— >^C, stars of spectral class C showing C/0> 1 (see e.g. Iben & Renzini 
1983; Smith & Lambert 1990). 

Another important consequence of TDU is the enrichment of the envelope in s-elements. 
The necessary neutrons for the s-process are released by two reactions: ^^C(q;, n)^^0, which 
provides the bulk of the neutron flux at low neutron densities (iV„ < 10^ cm~^), and 
^^Ne(a, n)^^Mg, which is activated at temperatures T > 3.0 x 10*^ K, providing a high peak 
neutron density (A^„ ~ 10^*^ cm^^) and is responsible of the production of s-nuclei controlled 
by reaction branchings (see Wallerstein et al. 1997; Busso, Gallino, & Wasserburg 1999 
and references therein). The abundances of s-nuclei are known to increase along the above 
mentioned spectral sequence, as the star gradually ascends the AGB. Evidence of this was 
provided during the last few decades by several studies on AGB stars of different spectral 
types: MS and S stars (C/0< 1) (Smith & Lambert 1985, 1986, 1990); SC stars (G/0~ 1) 
(LLoyd-Evans 1983; Abia & Wallerstein 1998, hereafter Paper I) and even post- AGB super- 
giants of spectral types A and F (Van Winckel & Reyniers 2000; Reddy, Bakker, & Hrivnak 
1999). 

The above studies found consistent enhancements of s-elements with respect to "nor- 
mal" red giants assumed as comparisons (or with respect to the Sun): from [s/Fe]~ -1-0.3 in 
MS giants to [s/Fe]> +0.5 in S stars"*^. The enrichment seems to continue along the AGB 
phase until the planetary nebula ejection. Indeed, post- AGB supergiants show high enhance- 
ments, [s/Fe]> +1.0. In this scenario, however, normal carbon stars (N-type)^ have still to 
find a place. This is so because of the complex spectra of N stars, which are so crowded 
with molecular and atomic absorption features (many of which unidentified) that abundance 
analysis has been strongly limited. 

The situation has not improved much with the advent of spectrum synthesis techniques. 
Indeed, the only abundance studies available to date are those by Kilston (1975) and Utsumi 
(1970, 1985), still based on abundance indexes and on low-resolution photographic spectra, 
respectively. These works suggested that N stars were s-element rich, showing overabun- 
dances by a factor of ten with respect to the Sun. As far as the AGB phase is concerned, 
this figure is not in direct contradiction with the accepted general scenario of s-nuclei en- 



^Wc; adopt the usual notation [X/Y]= log (X/Y)programstar- log (X/Y)comparisonstar for the stellar value 

of any clement ratio X/Y. 

^There are carbon stars of R-type and J-type, which probably owe their carbon enhancement to a mech- 
anism different than the TDU. These stars are not significantly enhanced in s-elements (Dominy 1985; Abia 
& Isern 2000) 



- 5 - 



hancement (see e.g. Busso et al. 1995). More likely, such high production factors might 
be difficult to reconcile with the enrichment in post-AGB supergiants showing similar, or 
even smaller, abundances. However, the large observational uncertainties for C-rich red gi- 
ants, and the lack of adequate model atmospheres has so far prevented a solid theoretical 
interpretation. 

In a previous work (Abia et al. 2001, hereafter Paper II) wc presented high-resolution 
spectroscopic observations for a sample of N stars focusing our attention on light s-elemcnts 
(Is), sited at the abundance peak near the neutron-magic number = 50, around the ^^Kr 
branching point of the s-proccss path. We showed how the analysis of the abundance ratios 
of Rb (a neutron- density sensitive element, see e.g. Beer & Macklin 1989) relative to its 
neighbors (Sr, Y and Zr) yields information on important details of the s-process mechanism 
operating, and on the initial stellar mass. We concluded that s-processing suggests low mass 
stars (LMS, M < 3Mq) as the hkely parents of C(N) giants. 

In LMS the major neutron source is ^^C, which burns radiatively in a tiny layer during 
the interpulse phase (Straniero et al. 1995) at relatively low temperatures (?» 8 keV), as 
a consequence of the formations of a ^^C-rich pocket in the intershell region. In the rarer 
AGB stars of intermediate mass (IMS) (M > 4Mo), ^^Ne would instead be favored as a 
neutron source, by the higher temperature in thermal pulses. Because of the very different 
neutron density provided by the two neutron-producing reactions, different compositions are 
expected from them, especially for the Is mixture. It is on this basis that in Paper II we 
drew our conclusions for the initial masses. 

In Paper II we also discussed the ^^C/^'^C ratios measured in the sample stars, showing 
that most of them cannot be explained by canonical stellar models on the AGB phase, 
requiring probably the operation of an ad-hoc mixing mechanism. This mechanism is often 
indicated with the term 'cool bottom process' (CBP). It is expected to occur in low-mass 
stars during the red giant branch and perhaps, also during the AGB phase (Wasserburg, 
Boothroyd, & Sackmann 1995; NoUett, Busso, & Wasserburg 2002). 

In the present work we have extended our study of s-element nucleosynthesis in N stars 
to the nuclei belonging to the second s-process peak: Ba, La, Ce, Nd, and Sm. We also 
re-analyzed the s-element abundances already derived in Paper II, by applying the spectral 
synthesis method. Together, these data are compared with recent models for s-processing 
in AGB stars at the metallicities relevant for our sample stars. Furthermore, using the 
infrared properties and Tc content, we discuss the possibility that our stars be extrinsic 
carbon stars, concluding that this is unlikely for several reasons. The structure of the paper 
is the following. In § 2 we present the characteristics of the stars and the spectroscopic 
observations. Section § 3 describes the method of analysis and the sources of error. Our 
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abundance results are reported in § 4 together with a comparison with similar studies in 
other AGB stars and with models of AGB nucleosynthesis. Finally, in § 5 we summarize the 
main conclusions that can be drawn from this work. 



2. Observations and analysis 

The spectra of our N star sample were collected at the Roque de los Muchachos (ORM) 
and Calar Alto (CAHA) observatories in several observational runs between the years 1998 
and 2001. In the former case wc used the 4.2 m William Herschel telescope and the 2.5 
m Nordic telescope with the UES and SOFIN spectrographs, respectively. The UES gives 
a maximum resolution of about 50000 between A = 450 and A = 900 nm, with some gaps 
between orders. SOFIN, working with the third camera configuration, reaches a maximum 
resolution of 180000 in 28 orders, covering about the same spectral range but with larger gaps 
between orders. At the CAHA observatory we used the 2.2 m telescope with the FOCES 
spectrograph which provides a resolution of about 35000 between A = 360 and A = 1000 
nm. In this case, an almost complete spectral coverage is obtained. Table 1 shows the 
sample stars and the epoch of the observations. Data reduction was performed following the 
standard procedures using the IRAF software package. Details on this are quoted in Abia & 
Isern (2000) and in Paper II. In the present work we add new N stars to the sample studied 
in Paper II, namely: RV Cyg, S Aur, SS Vir, T Dra, TT Cyg, U Hya, U Lyr, VY Uma, 
V429 Cyg and W CMa. 

2.1. Stellar parameters and luminosities 

The main parameters of the stars were derived in a way similar to the one of Paper 
II; they will not be discussed again here (see that paper for details). For all the stars we 
adopted a typical gravity of log g = 0.0 and a microturbulence of ^ = 2.2 km s^^. As in 
previous works, we used the grid of model atmospheres for carbon stars computed by the 
Uppsala group (see Eriksson et al. 1984 for details). The absolute abundances of CNO, the 
C/0 ratio (which mainly determines the global shape of the spectrum in carbon stars) and 
the ^^C/^^C ratios were derived as discussed in Paper II. Table 2 shows the final parameters 
adopted for each star. 

In order to estimate the absolute bolometric magnitudes Mboi of our stars (Table 1) 
we used the empirical relationships between Mboi and Mk for carbon stars obtained by 
Alkniss et al. (1998). The average K-magnitude values (the stars studied are variable!) 
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were taken from the Two Micron Sky Survey (Neugebauer & Leigthon 1969), available at 
the SIMBAD data base. Actually, though recent catalogs of near infrared images and point 
sources have become available (namely: 2MASS, Skrutskie et al. 1997, and DENIS, Epchtein 
et al. 1999), our sources are too bright for them, and are saturated. Recent observations 
at infrared wavelengths have been collected by the recent fifth edition of the Catalog of 
Infrared Observations, (Gezari et al. 1999); a few of our stars are present in the compilation, 
and the colors derived are compatible with those used here within the uncertainties. K 
magnitudes were corrected for interstellar extinction according to the Galactic extinction 
model by Arenou, Grenon, & Gomez (1992) with Ar/Av = 0.126 (CardeUi et al. 1989). 
Distances were calculated from the Hipparcos Intermediate Astrometric Data (van Leeuwen & 
Evans 1998). The Hipparcos parallaxes for carbon stars generally show much lower precision 
than those calculated for hotter main sequence stars. Parallaxes for carbon stars of different 
spectral types have been partially reprocessed (Knapp, Pourbaix, & Jorissen 2001; Pourbaix, 
Knapp, & Jorissen 2002; Bergeat et al. 2001, Bcrgeat, Knapik, & Rutily 2002a, b). However, 
the bolometric magnitudes showed in Table 1 have still to be taken with some caution. 
In any case, it is remarkable that the overwhelming majority of the stars have Mboi ^ 
—7.1, the maximum theoretical luminosity of an AGB star when its core mass reaches the 
Chandrasekhar mass limit, according to core-luminosity relation by Paczynski (1970). Our 
sample stars have in average Mboi ~ —5. but a considerable number have lower luminosity. 
This figure coincides with recent luminosity studies of carbon stars in the galaxy (Bergeat, 
Knapik, & Rutily 2002a) and in the nearby stellar systems (Westerlund et al. 1991, 1995) 



2.2. The list of lines 

Utsumi (1970) first noted that in N stars the blend effect of molecular bands is not 
very serious in the spectral regions A 4400 — 4500 and A ^ 4750 — 4950 A. He also found 
that many lines of s-process and rare-earth elements present in these spectral regions are of 
moderate intensity, and therefore they are useful for abundance analysis. Here, however, we 
have concentrated our analysis in the A = 4750 — 4950 A spectral range since the radiation 
fiux received from many carbon stars below 4500 A is very limited, because of the presence 
of strong C3 molecular absorption features below that wavelength limit. In Papers I and 
II we used the equivalent width versus curve of growth technique to derive heavy element 
abundances trough a careful identification and selection of the lines (see these papers for 
details). In the present study however, we decided to use the spectral synthesis technique 
because most of the lines from s-elements with high atomic mass (Ba, La, Ce, Nd and Sm, 
i.e. hs elements) present in this spectral range are heavily blended. This also allows us to 
verify the results of Paper II for the Is elements (Sr, Y, Zr), re-deriving their abundances 
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independently. This check confirmed the abundances deduced in Paper II (within the error 
bars, with the exception of a few cases, for reasons that are discussed below). This gives an 
important confirmation of the previous analysis technique. 

The spectral window studied, while being less affected by molecular blends than others, 
is nevertheless not free from such problems. Lines from the C2, CH and CN molecules have 
therefore to be considered. Transitions from the ^^€2, ^^C^''^C and ^^C2 (A-X), (B-A), (D-A) 
and (E-A) systems were found in the Kurucz database (cfaku5.harvard.edu/linelists.html). 
From the work of de Lavcrny & Gustafsson (1998) and Kurucz data, no C2 lines from the 
Phillips red system were predicted to be of importance in the studied spectral ranges. ^^C^^N, 
i3q14]\j-^ i2q15]\j- i3q15]\j- jj^es were included. The CN line lists were prepared along the 
same fines as the TiO fine list of Plez (1998), using data from Ito et al. (1988), Kotlar, 
Field, & Steinfeld (1980), Cerny et al. (1978), Rehfuss, Suh, & Miller (1992), Prasad & 
Bernath (1992), Prasad et al. (1992), Bauschlicher, Langhoff, & Taylor (1988) and Larsson, 
Siegbahn, & Agren (1983). Programs by Kotlar were used to compute wavenumbers of 
transitions in the red bands studied by Kotlar, Field, & Steinfeld (1980). Lines from the 
violet system were found to be very weak but were nevertheless considered. The absorption 
by the C2, ^^C^^N, ^^C^^N and ^^C^^N lines are found to slightly decrease the theoretical 
continuum by a few percent without affecting the line depths. In that spectral range, the 
main contributor from these two molecules is ^^C-'^^N. We adopted a solar ^^N/^^N in all the 
spectral synthesis. Several CH lines also play a role in the 4750 — 4950 A spectral range. The 
LIFBASE program of Luque & Crosley (1999) was used to compute these line positions and 
gf -values. Excitation energies and isotopic shifts (LIFBASE provides only fine positions for 
^^CH) were taken from the line list of J0rgensen et al. (1996). In this way, a very good fit 
of CH lines in the solar spectrum could be obtained, with the exception of very few lines 
which we removed from the list. Obviously, the line list also comprises atomic line data 
found in the VALD database (Kupka et al. 1998) and in Kurucz's data (CD-ROM No. 23) 
although, for the more intense atomic lines (mostly metals), we derived solar gf -values using 
the Holweger & Miiller (1974) model atmosphere and solar abundances by Anders & Grevesse 
(1989). The gf -values derived in such a way showed an excellent agreement with the solar 
values derived by Thevenin (1989, 1990). For a sample of about 100 lines in common we 
found a mean difference of 0.15 ±0.2 dex. Our final line list includes more than 20000 atomic 
and molecular lines. Despite of this, only a few lines were selected as indicators of the hs 
element abundances in our stars. These lines are shown in Table 3. The metallicity [M/H] 
of our sample stars was derived from selected metallic lines, "M" indicating an average value 
of Ca, Ti, V, Fe and Co abundances (see Table 3 in Paper II). Using this average changes 
only slightly the values [Is/Fe] and [hs/Fe] derived. 
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2.3. Data analysis, peculiar cases, comparison with previous work 

The hs clement abundances were derived from spectrum synthesis fits to the observed 
spectra. The theoretical spectra were computed in LTE and the abundance analysis was done 
line-by-line relative to the star WZ Cas (C9,2J;SC7/10) which was used, as in Paper II, as the 
reference star. WZ Cas is a carbon star with a C/0 ratio very close to unity with near solar 
metallicity and with very small or no s-element enhancement: [h/Fe]= 0.02 ± 0.20, where h 
(heavies) includes Sr, Y, Zr, Nb, Ba, La, Ce and Pr (see Lambert et al. 1986; Abia & Isern 
2000). This comparative analysis might reduce systematic errors due to wrong 51/- values and 
possible departures from LTE (Tomkin & Lambert 1983). All the abundance ratios derived 
(see below) are relative to this star. The final abundance derived for each element is the 
mean of the abundances derived from each line, except for Ba whose abundance is derived 
from only one line. In deriving these means, we consider twice the abundances derived from 
particular lines which we believe are better abundance indicators for a given element because 
of their larger sensitivity to abundance variations and/or because they are placed in a less 
crowded region. These fines are: the A4809 Lall, A4773 Cell and A4811 Ndll lines (see Table 
3). 

Figure 1 shows an example of a theoretical fit to the observed spectrum of the star Z Psc 
in the region between 4800 — 4830 A, where most of the available s-element fines are located. 
Figure 2 is a similar example in the region around the 4934 A Ball line. It can be seen in 
both figures that the observed spectra are generally well fitted by spectral synthesis, despite 
the many features still not accounted for. This last fact is a clear evidence of the remaining 
incompleteness of the line list used, despite the several thousands of lines included. 

As noted above, in the present work we have re-derived with the spectral synthesis 
technique the abundances of Sr, Y and Zr obtained in Paper II; this was permitted by the 
fact that most of the available lines for these elements are located in the same spectral range 
4750-4950 A studied here (see Table 3 in Paper II). We find a good agreement (within typical 
error bars) with the values found for these elements in Paper II; hence, our conclusions are 
not modified by the present analysis. Only for one star, V460 Cyg, the abundances of Sr, Y 
and Zr are significantly modified. The reason of this lays in a bad interpolation done in Paper 
II over the grid of model atmosphere parameters for this star. There is also a modification 
of the Rb and Zr abundances in Z Psc with respect to Paper II. The Rb abundance shown 
there (Table 2) is just a misprint, the correct value being [Rb/M]= 0.6 instead of 0.9. The 
revised [Zr/M] ratio here is +1.0 dex, which is, nevertheless, in agreement within error bars 
with the previous value in Paper II (+0.8 dex). In the present work we have also been able 
to estimate the Sr abundance in this star using the Sr I line at 4811.88 A: we recall that the 
Sr abundance in Z Psc was not derived in Paper II. 
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in Paper II we also discussed the increased difficulties found in the analysis of certain 
N stars, in particular lY Hya, VX Gem and V CrB (see that paper for details). In fact, 
the abundance results in such stars were considered with a lower weight in comparing with 
theoretical models. Unfortunately, here the list of difficult stars has to be extended. Indeed, 
T Dra, U Lyr (both Mira variables), S Aur and V429 Cyg show spectra so smeared out 
that we had no choice but decide to give up their analysis. For SS Vir and V Oph (also 
Mira variable), the spectra are also smeared out, and in fact they provide an insufficient 
agreement between observed and theoretical spectra. In particular, all Miras in our sample 
show difficult spectra. It is not clear why spectral lines in the Mira variables look in this 
way. Mira variables have generally thick circumstellar dust shells, whose contamination due 
to thermal emission by dust might be responsible for this problem. The gas motions due to 
pulsation might contribute to the broadening of spectral lines. Finally, the star V781 Sgr 
shows a weak central emission at the core of the Na D lines. This is usually interpreted as 
due to the presence of a chromosphere. The influence of a chromosphere in the derivation of 
abundances in carbon stars is uncertain and is beyond the scope of this work; therefore, we 
have included also this star in the group of difficult N stars. Summing up, the abundances 
derived for lY Hya, VX Gem, V CrB, SS Vir, V Oph and V781 Sgr have to be considered 
with considerable caution. To emphasize this, in the figures below, we have plotted with 
smaller symbols the abundance ratios derived for them. 

The stellar atmosphere parameters (Teff, gravity, etc) for N stars are difficult to de- 
termine with accuracy (Bergeat, Knapik, & Rutily 2001, 2002a). Note for instance that 
they are variable stars, therefore the effective temperature changes during the pulsation. 
The microturbulence in variable stars depends on the wavelength in the spectrum (due to 
differences in depth of formation), and changes with the pulsational phase. Furthermore, 
the existence of strong waves, shock fronts, and inhomogeneities is the rule rather than the 
exception in their atmosphere (sec Gustafsson & J0rgenscn 1994 for discussion). In Table 4 
we show the sensitivity of the abundances derived to typical uncertainties in the atmosphere 
parameters. Since most of the lines used in the present work correspond to singly ionized 
elements (see Table 3), an important source of uncertainty is the error in the gravity. The 
abundances derived are rather sensitive to the microturbulence parameter, indicating some 
saturation of the lines used in the analysis. However, the hs element abundances arc mostly 
insensitive to uncertainties in metallicity, CNO absolute abundances and the ^^G/^^C ratio. 
Certainly, deviations from LTE might play a signiflcant role, although they are reduced in 
a differential analysis as is done here. If we consider the uncertainty due to the placement 
of the spectral continuum (~ 5%), adding quadratically all the sources (non systematic) of 
error, we found a typical uncertainty in the absolute abundances of about ±0.2 dex for Ba, 
±0.4 dex for La, Nd, and Sm and ±0.45 dex for Ce. 
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Furthermore, we have to consider the random error in the abundances of elements 
represented only by few lines, as in our case. Wc can estimate this random error only 
in stars where we derived individual element abundances from more than three lines. In 
this case we found a dispersion in the range 0.1 — 0.3 dex. Nevertheless, when computing 
the [hs/M] ratios, errors are lower because some sources of error affect in the same way 
the metallicity (see Table 4 in Paper II). Adding up all these contributions, we estimate a 
typical total error (non systematic) in the [hs/M] ratios ranging from ±0.3 (Ba) to ±0.45 
dex (Ce). A similar discussion can be done concerning the uncertainties in the derivation 
of the Is element abundances and metallicity (see Paper II for details). Considering this, 
we estimate a typical error of about ±0.3 dex in the [hs/ls] ratios. In Table 5 we show 
the abundance ratios with respect to the metallicity derived in our stars, referred to the 
comparison star WZ Cas. The abundances in Table 5 represent the first detailed analysis of 
the s-clement composition in N stars. In this table we show the previous results (revised) 
obtained in Paper II for the Is elements, together with the new values derived here for the 
hs elements, including some additional stars (see §2). We confirm our previous finding that 
N stars are of near solar metallicity: <[M/H]>= —0.02 ± 0.10. Only two stars seem to be 
of low metal content (lY Hya and V CrB); however, we remind the difficulties found in the 
analysis of these objects (see above). Considering the abundances of Y and Zr as the only 
good tracers of the Is element enhancement^, our objects have <[ls/M]>= +0.67 ± 0.24. In 
the same way, considering Ba, La and Nd as representative of the hs element enhancement 
(Ce and Sm are derived only in a few stars; in many of them we set just upper limits), the 
mean is <[hs/M] >= +0.52±0.29. This figure represents moderate, but significant s-element 
enhancement, which is clearly smaller than that obtained previously by Utsumi (1985). We 
refer the reader to Paper II for the detailed discussion made there about the explanation of 
this difference. 



3. Tc and possible binarity 

The radioactive nucleus ^^Tc (ri/2 ~ 2 x 10^ yr) is produced along the s-process path. 
Because the AGB phase spans only a few half-lives of this element, once dredged- up into the 
envelope, it should remain present in a certain percentage throughout the AGB lifetime. 

Therefore, the detection of Tc in the atmosphere of an AGB star (Merill 1952) is an 
undeniable evidence of the s-process working in situ {intrinsic AGB). When, on the contrary. 



■^Because of the difficulties in deriving Sr abundances, we discard this element when computing the mean 
low-mass s-element abundances Is (see discussion in Paper II). 
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Tc presence can be excluded, the star most probably owes its carbon enrichment to a mass 
transfer episode from an AGB companion, extrinsic AGB. Unfortunately, the detection of Tc 
in AGB stars is a difficult task. The most interesting Tc lines are located around A ~ 4260 
A, where AGB stars are not very bright. 

Indeed, in most N stars this spectral region is not accessible to observations because 
of the strong flux depression at these wavelengths. There arc. however, a couple of inter- 
combination Tc lines that may be used for the analysis. One of them is the Tel line at 
A 5924.47 A, which is however weak and not free from blends. We have extended the 
analysis of Tc in N stars made in Paper 11 using this line in the additional stars included in 
this work. We have followed the same technique as described in Abia & Wallcrstcin (1998). 
In Table 5 we show the results of this search. As in Paper II, a "yes" entry in this table 
means that the best fit to the 5924 A blend is obtained with a non-zero Tc abundance. A 
"no" entry means that a synthetic spectrum with no Tc does not significantly differ from 
another one with a small Tc abundance and, finally, a "doubtftil" entry means that the 
5924 A blend can be fitted better with a small abundance of Tc, but the Tc absorption in 
the blend is not so clearly appreciated as in the case of the Tc "yes" stars. If we consider 
the "doubtful" stars as yielding a detection, which is probably the case, Table 5 shows that 
about 62% of our targets are Tc-yes N stars and thus, can be considered as bona fide intrinsic 
carbon stars. However, because of the difficulty in the analysis of this Tc line (weak, and 
in a very crowded zone), a non-detection does not necessarily mean that Tc is absent: we 
do not have, unfortunately, a clear-cut possibility of excluding it. Probably, a study of the 
resonance lines at A ~ 4260 A, when feasible, could yield more secure conclusions. Thus, 
additional tests are necessary to establish the extrinsic or intrinsic nature for the rest of the 
sample. 

A further help can luckily come from an analysis of infrared colors. This idea (Jura 
1986; Little- Marenin & Little 1988) arises from the fact that intrinsic AGB stars (in our case 
N stars) are high mass losers and show excess infrared radiation from circumstellar dust, 
whereas extrinsic N stars would be first-ascent red giants with mass loss rates smaller by two 
orders of magnitudes at least. This was already used by Jorissen et al. (1993) to distinguish 
between extrinsic and intrinsic S stars. These authors showed that the overwhelming major- 
ity of Tc-deficient S stars have a ratio R = F {12 jim) / F {2.2 jim) < 0.1, where F(12/xm) and 
F(2.2/im) are the 12 //m and 2.2 //m fluxes from the IRAS Point Source Catalog and from 
the Two-Micron Sky Survey, respectively. On the contrary, most Tc-rich S stars (intrinsic) 
have R > 0.1. We have performed a similar study in our sample of carbon stars and the 
result can be seen in the histogram of Figure 3. We have considered in this plot the doubtful 
stars as showing Tc in their spectra, because of the reasons commented above. From this 
figure it is obvious that most of our stars have infrared excess R> 0.1 and only a few show 



-13- 



lower values (in fact in the very narrow range ~ 0.08 — 0.1). It is clear that, contrary to what 
happens for the S stars (sec Jorissen et al. 1993), there is no distinction in the distribution 
between Tc-yes and Tc-no carbon stars. We might then conclude that probably all carbon 
stars in our sample are in fact of intrinsic nature and the quoted Tc-no carbon stars are just 
a consequence of the difficult and uncertain analysis of the Tc I 5924 A blend on which this 
study is based^. 

Similarly to what happens in S stars, a difference between Tc-yes and Tc-no carbon stars 
might be also apparent in their temperature class. The large uncertainty in the derivation of 
the effective temperature in carbon stars docs not allow us to address this question directly 
from the estimated T^g values. The point is that Tc-yes carbon stars should be cooler and 
more evolved and hence should (on average) have more circumstellar matter to account for 
the infrared excess. This question can be addressed by studying the K— [12] vs. K— [25] 
color-color diagram. The color index K—[i] is defined as K— 2.5 log(620/F(i)), F{i) being 
the (non color corrected) flux in the i band taken from the IRAS Point Catalog Source. 
Again, we do not find a distinction between Tc-yes and Tc-no stars (see Figure 4). All the 
stars are located in the region of black bodies with T < 3000 K. These temperatures are 
indeed expected in evolved stars during the AGB phase of evolution (note than some stars 
have black body temperatures lower than 2000 K). Higher black body temperatures would 
instead be expected in extrinsic (RGB) carbon stars. 

A better discrimination will be possible in this field, when extensive colors for C-rich 
and 0-rich AGB stars from modern infrared cameras will become available, providing data 
with a higher precision than possible from IRAS (see preliminary results by Marengo et al. 
1999; Busso et al. 2001b). 

Finally, we can also address the problem of the nature of our sample stars (i.e. those 
being extrinsic or intrinsic AGBs) through a comparison with theoretical models of element 
enrichment in the two classes of objects, to look for parameters helping in the discrimination. 
There are several questions for which we can, in this way, find useful answers. Is it possible, 
at nearly solar metallicity, to produce an extrinsic carbon star (C/0> 1) from the transfer of 
carbon-rich material in a binary system (most probably through the wind accretion mecha- 
nism, Jorissen &; Mayor 1988)? What should be the C/0 ratio in the transferred material so 
that the mass-receiving companion remains C-rich? How much these requirements depend 
on the stellar mass and on the metallicity of the accreting star and of the donor? We shall 



''Note that, as mentioned in Paper II, in some stars we were able to study simultaneously the stronger 
resonance Tcl lines at A ~ 4260 A. In all the stars where this was possible, we confirmed the detection found 
from the 5924 A blend. 
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come again on these problems in the next section, where we present a comparison of our 
observations with the results of nucleosynthesis and mixing models for intrinsic and extrinsic 
AGB stars. 



4. Comparison with models of AGB nucleosynthesis 

In order to understand the way in which the observed abundances were created and 
subsequently mixed into the stellar convective envelopes, we made use of the same stellar and 
nucleosynthesis models already described in Paper II. They start from full stellar evolution 
calculations, computed with the FRANEC evolutionary code, spanning the metallicity range 
from solar to 1/20 solar, and the mass range from 1.5 to 7 Mq (Straniero et al. 1997, 2000). 

For our purposes, we adopt the 1.5 Mq model as representative of LMS, since in Paper 
II we already concluded that most N stars are of low mass (M < 3Mq). Notice that the 
theoretical expectations for a 3 Mq stellar model mimic quite closely the 1.5 Mq model: 
differences in s-element logarithmic abundances arc found to be below 0.1 dcx. We have 
computed models without mass loss, from which an asymptotic value of the parameter 
A = AM(TDU)/AMh has been derived, namely 0.24 for the 1.5 Mq model. On this basis 
a post-processing has been performed to derive the detailed nucleosynthesis products in 
the He intershell and the chemical modifications of the envelope (Gallino et al. 1998). The 
duration of the thermally pulsing AGB phase is mainly determined by the assumed mass loss 
rate. Hence, mass loss was included in the nucleosynthesis calculations trough the Reimers's 
(1975) phenomenological law, choosing the value 0.7 for the free parameter 77. In each 
interpulse/pulse cycle we compute the detailed nucleosynthesis pattern that is established 
in the intershell zone, induced first by the neutron release from radiative ^^C burning in the 
interpulse phase, and then by the activation of the ^^Ne neutron source in the convective 
instability of the He shell. One cycle after the other, we follow the penetration of envelope 
convection through TDU that mixes part of the processed material to the envelope; the whole 
sequence is repeated till our models experience the end of TDU. The procedure assumes 
that some form of proton penetration below the convective border can occur at dredge up, 
polluting the radiative He- and C-rich layers, so that at the restart of H-shcll burning the 
main neutron source ^^C can be formed. For the problems and the uncertainties involved in 
this assumption, as well as for the techniques used in the computations of neutron captures 
and the nuclear parameters used, see Paper II and Gallino et al. (1998). There, in particular, 
we described the way in which our model sequences are created, and the efficiency of ^^C 
burning assumed. As in Paper II, we shall refer with the term ST (standard) to models in 
which the mass of ^^C burned per cycle is 4xl0~^ Mq in LMS. Then this reference abundance 
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is scaled upward and downward to span a wide range of possibilities in the nucleosynthesis 
efficiency. 

Figures 5 to 7 show a synthesis of our predictions: they present, as a function of the 
initial metallicity [M/H], the average photospheric abundance ratios [hs/M], [Is/M], [hs/ls] 
at C/0 —1, in our 1.5 Mq model, assumed as representative of the LMS class. 

Inspection of Figures 5, 6 and 7 provides a quite satisfactory agreement between obser- 
vations and theoretical models, inside the wide spread of neutron release efficiency considered 
and taking into account the observational uncertainties of N stars. This is in agreement with 
previous findings on the average s- process efficiency by Busso et al. (2001a), derived by 
other classes of AGB stars and their descendants at various metallicities. 

A rather comprehensive picture of the s-enrichment in N stars is provided by plots in 
which observations and theory are compared using the s-element enhancement (indicated by 

e.g. Is) and a parameter sensitive to the neutron exposure (eg [hs/ls], as shown in Busso 
et al. 1995, 2001a). This kind of comparison is performed in Figure 8, where N stars are 
shown together with other known classes of intrinsic and extrinsic AGB stars. As the figure 
shows, the spread in the s-proccss efficiency (as monitored by the N("'^^C)/N(^^Fe) ratio in 
the pocket), is in this way constrained better, being rather limited. Another result is that N 
stars are very close to S stars for what concerns the efficiency of mixing into the envelope. In 
the models, it is indeed often sufficient to make one extra pulse to change the composition 
of an S star into that of a N giant, establishing a C/0 ratio larger than unity. This has large 
consequences on the appearance of the atmosphere and of the circumstellar envelope (chem- 
istry, colors, type of dust formed all change considerably), but the corresponding change 
in the s-process enhancement may be rather limited. Since post-AGB stars usually show 
larger s-elemcnt enhancements than those derived here in N-stars (see § 1), and because our 
sample stars show a C/0 ratio very close to one, we might speculate the stars analyzed here 
have just become carbon stars. As the star evolves in the AGB phase and TPs and TDUs 
continue operating, more carbon and s-elements are dredged-up into the envelope. Eventu- 
ally the C/0 ratio exceeds unity by a such large amount that the star becomes obscure at 
optical wavelengths due to the formation of dust. This scenario might explain the apparent 
gap in the level of s-element enhancement and C/0 ratio observed between C(N) stars and 
post-AGB stars: we cannot "see" carbon stars with a C/0 ratio largely exceeding unity. We 
have to outline once again that our conclusions are for the moment somewhat hampered by 
the still large observational uncertainties for N stars, whereas typical error bars for the other 
classes of s-enriched stars in the Galactic disc are of in general lower (0.20 to 0.25 dex, see 
e.g. Busso et al. 2001a). 

The interpretation of the N star abundances presented so far in the framework of LMS 
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evolution along the AGB is confirmed, and even strengthened, by an analysis of the detailed 
element distributions in the observed stars, from the data of Table 5. In order to show 
this we have plotted in Figure 9 a sample of fits to individual stars, obtained by comparing 
their observed composition with that yielded by the models in the photosphere, at C/0 =1. 
The plots show a complete agreement, and wc believe this provides a full confirmation that 
most, if not all, N stars of the Galactic disk are indeed of low initial mass. In fact, from 
the discussion presented in Paper II we know a priori that plots made similarly to Figure 9, 
but using models from IMS, would be often incapable of reproducing at least the "Is" zone, 
yielding in general a too high Rb/Sr ratio. This conclusion for the mass of carbon stars is 
in line with the present independent understanding of the different classes of carbon stars, 
provided by astrometric, kinematic, and photometric criteria (Bergeat, Knapik, & Rutily 
2002a, b). 



4.1. Can any of our sources be an extrinsic C-star? 

Using the models outlined in the previous subsection, we can add something to the 
discussion of the nature of our sources. This can be done on an independent basis with 
respect to the (already analyzed) presence of Tc and infrared properties. Having in mind 
the model compositions discussed above, we can ask ourselves which predictions for s-element 
and carbon abundances would be derived in case some of our sample stars were extrinsic 
AGBs. In order to study this possibility, we have first to consider that the resulting C/0 
ratio after mass transfer is constrained by the observed C/0 value, which is very close to 
1. Then we must take into account the original C and O content in the companion. This 
depends on metallicity, and at the moment of mass transfer it has probably not changed from 
the original CNO composition. Indeed, the long duration of core H burning (that precedes 
any significant abundance change at the surface) makes this phase the most likely for the 
occurrence of mass accretion. Subsequently, partial ^^C dilution is produced by the first 
dredge up (FDU), during the RGB phase; this also increases the abundance of ^^C (^^O is 
only slightly reduced by FDU). 

For a star reaching the RGB with a solar mixture of CNO isotopes, FDU would change 
CN abundances by a significant fraction (e.g. -"^^C would decrease by 13% to 40% in solar 
metallicity models spanning the range 1-5 Mq). Its impact is however marginal when the 
star has a high carbon enhancement at the surface. Furthermore, in case of a lower-than- 
solar metallicity, one must take into account that the initial carbon abundance decreases 
with Fe, whereas [0/Fe] shows a concomitant enhancement. One can adopt general rules for 
this enhancement, deducing it from observations of large databases of dwarf stars at different 
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metallicities, e.g. as described in Israelian et al. (2001). Note that while [0/Fe] vs. [Fe/H] 
still remains controversial at [Fe/H]< —1, this does not affect our discussion, because the 
typical metallicity of our sample stars is much higher. 

We then modeled mass transfer by diluting the AGB envelope composition, allowing 
for different "mixing histories" , which take into account the different possibilities described 
above. In all our attempts wc found that, for metallicities not far from solar, it is impossible 
to produce an extrinsic carbon star. This gives a strong support to the indications already 
obtained from the observational criteria. 

To make this result easily understood, we can express the mass transferred by means 
of a semi-analytic formula. If mass transfer occurs while the secondary component is on 
the main sequence, the material accumulates onto a thin radiative region, whose chemical 
composition becomes that of the accreted material from the AGB companion and is therefore 
characterized by a C/0 ratio in excess of unity. At FDU, a homogeneous composition is 
established by mixing the accreted material over more internal layers, where the original 
chemical composition, partly modified by H burning, has been preserved. Then, we can 
define a dilution factor / as the ratio of the transferred mass (M^q^) to the fraction of the 
envelope having the original composition (M^^™^). By imposing that, in the photosphere 
of the secondary star, the C/0 ratio remains equal to 1 even after FDU, one can see that: 

/ = (1 - 0.28/10-°-^[^^/Hl)((C/O)JfG*™u _ 

Here we have assumed that mass transfer occurs at the last TDU episode, in order to 
allow the maximum possible carbon enhancement of the secondary component, thus favoring 
the formation of an extrinsic carbon star. According to the FRANEC models, TDU ceases 
when the envelope mass drops below a critical limit, namely ~0.5 M©; from then on, while 
mass loss can further decrease the envelope mass, its C/0 ratio remains frozen. In Table 6 
we report, for the 1.5 Mq model at different metallicities (column 1), the final C/0 ratio 
of the intrinsic AGB star (primary component) at the last thermal pulse that experienced 
TDU ((C/0)J^G™^) ^-iid the corresponding dilution ratio /, as derived by means of the 
formula reported above. Obviously, there is a maximum limit in the transferred mass, which 
corresponds to the accretion of the whole residual mass of the AGB envelope (i.e. ~0.5 
Mq). In order to give an estimate for such a limit, we further assumed that the mass 
of the companion star, after the mass accretion, is the lowest compatible with the age of 
the Galactic disk, let us say ~1 Mq. Again, this is intended to favor carbon enrichment, 
as it roughly corresponds to the minimum mass of the convective envelope at the FDU. 
Correspondingly, in the most favorable conditions, we obtain an upper limit for the dilution 
factor / of 1.5. We emphasize that this limit corresponds to a very extreme case and we 
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expect that a substantially lower value of the / parameter should be more reahstic. For 

example, if the mass of the secondary component after the mass accretion episode were 1.5 
Mq (corresponding to an initial mass of about 1 Mq), the limiting value of the / parameter 
would become 0.7. In any case, from a comparison with the calculated / values reported 
in Table 6, it is clear that it would be practically impossible to form extrinsic carbon stars 
for solar metallicities. We suggest that the maximum metallicity of an extrinsic carbon star 
should be of the order of [Fe/H] 0.4 to -0.3. 

Thus, even considering the uncertainty in the metallicity of the stars considered in the 
present work (±0.3 dex), this simple result prevents the vast majority of them from being 
considered as extrinsic carbon stars. 

Note that if a carbon star is formed as a consequence of the mass transfer process, as is 
indeed the case at lower metallicities, then the behavior of [Is/Fe] versus [Fe/H] would not 
be distinguishable from the trend seen for intrinsic AGBs. This is related to the fact that 
mass transfer and the following FDU behave much hke a form of dilution, as in the case of 
TDU, and we have no way to distinguish which of the two possible forms of dilution the 
s-enriched material has undergone (see Busso et al. 2001a, and in particular Figures 11 and 
12 there, for a discussion of this in the case of Ba and CH stars). 

In conclusion, both the observational and theoretical analysis lead us to conclude that 
the objects studied here are most probably intrinsic carbon stars, presently evolving along 
the TP-AGB. Exceptions may exists, but this should be tested with more direct methods, 
like radial velocity variations and their effects on the spectra. At the moment, for none of 
the N stars studied here there is any evidence of such variations. 

5. Conclusions 

In this paper we have presented for the first time a rather complete sample of s-process 
abundance measurements in N-type carbon stars, extended to light and heavy species, across 
the s- process peaks with neutron numbers = 50 and N = 82. Making use of a large 
line list, and of up-to-date model atmospheres for carbon rich giants, we have deduced 
the s-element abundances through spectrum synthesis technique. In doing so, we have 
also verified previous results (Paper II) for the light s-elements across the *^"''Kr branching 
point of the s-path. N stars turn out to be characterized by s-element abundances very 
close to, or slightly higher than, those found for S stars. Compared with the only analysis 
performed before (Utsumi 1985), our abundances are clearly smaller; they are smaller also 
with respect to those of post-AGB supergiants (Reddy, Bakker, & Hrivnak 1999; Van Winckel 
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& Reyniers 2000), and this gives rise to a general picture in which the surface enrichment 
continuously increases along the evolutionary sequence producing MS, S, and N stars, and 
subsequently yellow post-AGB supergiants (Rcddy et al. 2002). We have compared our data 
with model envelope compositions obtained from previously published calculations of AGB 
nucleosynthesis and mixing. Good agreement is obtained between low mass star models 
and s-element observations. Several pieces of evidence (from the detection of Tc lines and 
infrared colors, to the theoretical relations between the initial metallicity, the mass, and 
the final C/0 ratio) lead us to conclude that most (or perhaps all) our sample stars are 
intrinsic TP- AGB stars, so that their abundances are locally produced by the occurrence of 
third dredge-up during the TP-AGB phases, and not generated by mass transfer in binary 
systems. 
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Table 1. THE PROGRAM STARS 



Star Spectral^ Var.^ Period^ Mboi^ Epoch Observatory^ 
type type (days) (JD+2450000) 



AO And 


C5,4 


SR 


346 


-5.2 


791 


CAHA 


AW Cyg 


C3,6 


SRb 


220 


-5.7 


947 


ORM 












1745 


CAHA 


EL Aur 


C5,4 


Lb 




-2.5 


791 


ORM 


HK Lyr 


C6,4 


Lb 




-6.0 


947 


ORM 












1015 


CAHA 


IRC -10397 


N 








1015 


CAHA 












1745 


CAHA 


lY Hya 


N 


Lb 




-5.2 


947 


ORM 


LQ Cyg 


C4,5 


Lb 






1745 


CAHA 


RV Cyg 


C6,4 

* 


SRb 


263 


-7.1 


1745 


CAHA 


RX Set 


C5,2 


Lb 






1015 


CAHA 


S Aur 


C5,4 


SR 


590 




745 


CAHA 


S Set 


C6,4 


SRb 


148 


-4.6 


1015 


CAHA 


SS Vir 


C6,3 


SRa 


364 




1745 


CAHA 


SY Per 


C6,4 


SRa 


476 


-3.1 


791 


CAHA 


SZ Sgr 


C7,3 


SRb 


100 


-2.1 


1015 


CAHA 


T Dra 


C6,2 


Mira 


421 


-2.0 


791 


CAHA 


TT Cyg 


C5,4 


SRb 


118 


-4.6 


1015 


CAHA 












1745 


CAHA 


TY Oph 


C5,5 


Lb 






1015 


CAHA 


U Cam 


C3,9 


SRb 


400 


-4.8 


791 


ORM 


U Hya 


C6,3 


SRb 


450 


-4.0 


947 


ORM 












1919 


ORM 


U Lyr 


C4,5 


Mira 


457 


-5.5 


791 


CAHA 


UU Aur 


C5,3 


SRb 


235 


-6.0 


947 


ORM 


UV Aql 


C6,2 


SRa 


385 


-7.0 


947 


ORM 












1745 


CAHA 


UX Dra 


C7,3 


SRa 


168 


-5.6 


947 


ORM 


V Aql 


C6,4 


SRb 


353 


-5.0 


1745 


CAHA 
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Table 1 — Continued 



Star 


SpcctraF 


Var.i 


Period^ 


Mboi' 


Epoch 


Observatory^ 




type 


type 


(days) 




(JD+2450000) 




V CrB 




Mira 


357 


-7.1 


1U15 


CAHA 












947 


(JRM 


V (Jpn 


Co, 2 


Mira 


298 


-5.9 


1015 


A TU A 

CAHA 












947 


(jRiVl 


VX Gem 


C7,2 


Mira 


379 


-3.6 


745 


CAHA 












1919 


URM 


VY Uma 


C6,3 


Lb 




-4.3 


1919 


CRM 


V429 Cyg 


C5,4 


SRa 


164 




745 


CAHA 












1745 


CAHA 


V460 Cyg 


C6,3 


SRb 


180 


-5.8 


1015 


CAHA 












1745 


CAHA 


V781 Sgr 


N 


Lb 






1015 


CAHA 


W CMa 


C6,3 


Lb 




-7.3 


1919 


ORM 


W Ori 


C5,4 


SRb 


186 


-5.6 


791 


CAHA 


WZ Cas 


C9,2J 


SRa 


186 


-6.2 


745 


CAHA 












791 


ORM 












1919 


ORM 


Z Psc 


C7,2 


SRb 


144 


-4.1 


745 


ORM 



^Data taken from the SIMBAD database. 

^Bolometric magnitudes are estimated from mean K magnitudes and HIP- 
PARCOS parallaxes (see text) . 

^CAHA: Centro Astronomico Hispano-Aleman de Calar Alto; ORM: Ob- 
servatorio del Roque de los Muchachos. 
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Table 2. DATA FOR PROGRAM STARS 



Star Teff Ref.^ C/0 ^'^C/^^C Ref.^ K-[12]^ K-[25]= 



AO And 


2970 


2 


1.02 


30 


1 


-1.80 


-3.16 


AW Cyg 


2760 


2 


1.03 


21 


2 


-1.59 


-2.92 


EL Aur 


3000 


6 


1.07 


50 


2 


-1.64 


-2.83 


HK Lyr 


2866 


2 


1.02 


10 


2 


-1.92 


-3.18 


IRC -10397 


2600 


6 


1.01 


20 


2 






lY Hya 


2500 


6 


1.02 


15 


2 


5.61 




LQ Cyg 


2620 


2 


1.10 


40 


2 


-1.98 


-3.19 


RV Cyg 


2600 


1 


1.20 


74 


3 


-2.87 


-1.59 


RX Set 


3250 


2 


1.04 


60 


2 






S Set 


2895 


1 


1.07 


45 


3 


-1.92 


-3.32 


SS Vir 


2560 


7 


1.05 


10 


2 


-2.60 


-1.34 


SY Per 


3070 


3 


1.02 


43 


2 


-1.69 


-2.86 


SZ Sgr 


2480 


2 


1.03 


8 


2 


-1.74 


-2.43 


TT Cyg 


2825 


7 


1.04 


30 


2 


-3.50 


-2.05 


TY Oph 


2780 


2 


1.05 


45 


2 


-1.89 


-2.95 


U Cam 


2670 


2 


1.02 


40 


2 


-1.44 


-2.57 


U Hya 


2825 


1 


1.05 


35 


2 


-3.00 


-1.86 


UU Aur 


2825 


1 


1.06 


50 


3 


-1.80 


-3.06 


TIV Anl 


2750 


2 


1 00^ 


1 9 


2 


-1 74 


-3 06 


UX Dra 


2900 


1 


1.05 


26 


2 


-2.30 


-3.58 


V Aql 


2610 


1 


1.16 


90 


2 


-1.74 


-3.17 


V CrB 


2250 


4 


1.10 


10 


4 


0.02 


-1.24 


V Oph 


2880 


2 


1.05 


11 


1 


-1.64 


-2.96 


VX Gem 


2500 


5 


1.01 


60 


2 


-1.46 


-2.93 


VY Uma 


2855 


1 


1.08 


44 


3 


-3.56 


-2.13 


V460 Cyg 


2845 


1 


1.06 


61 


3 


-2.05 


-3.43 


V781 Sgr 


3160 


2 


1.02 


35 


3 


-3.22 


-1.86 


W CMa 


2880 


1 


1.09 


53 


3 


-2.92 


-1.98 


W Ori 


2680 


1 


1.20 


79 


3 


-3.04 


-1.69 


WZ Gas 


3140 


3 


1.005 


4 


3 


-3.66 


-2.39 


Z Pse 


2870 


1 


1.01 


55 


3 


-3.55 


-2.39 
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Table 2 — Continued 



Star Teff Ref.^ C/0 ^^C/^^G Ref.*^ K-[12]'^ K-[25]^ 



^References for T^^: (1) Lambert ct al. (1986); (2) Ohnaka & Tsuji (1996); 
(3) Dyck, van Belle, & Benson (1996); (4) Kipper (1998); (5) Groenewegen et 
al. (1998); (6) For these stars we have estimated the effective temperature by 
comparing their spectra with the temperature sequence spectral atlas for C-stars 
created by Barnbaum et al. (1996); (7) Bergeat, Knapik, & Rutily (2001) 

•^Source for the ^^C/i^C ratio: (1) Ohnaka & Tsuji (1996); (2) Derived in this 
work or Paper II; (3) Lambert et al. (1986); (4) Kipper (1998). 

^The 12 and 25 //m fluxes are taken from the second edition of the IRAS PSC 
anf the K magnitudes from the Two- Micron Sky Survey indicated in Claussen et 
al. (1987). 
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Table 3. SPECTROSCOPIC DATA 



Element A(A) x(eV) log gf Reference ^ 



Ba II 


4934, 


.10 


0. 


,00 


-0.150^ 


1 


La II 


4804, 


.04 


0. 


,23 


-1.530 


5 


La II 


4809, 


.02 


0. 


,23 


-1.267 


2 


La II 


4748, 


.73 


0. 


.93 


-0.600 


2 


Ce II 


4773, 


.96 


0. 


.92 


0.010 


3 


Ce I 


4820, 


.00 


0. 


,27 


-0.559 


4 


Ce I 


4820, 


.61 


0. 


,17 


-0.866 


4 


Ce I 


4822, 


.54 


0, 


,00 


-0.420 


4 


Nd II 


4811, 


.35 


0. 


,06 


-0.570 


s 


Nd II 


4817, 


.18 


0. 


,74 


-0.760 


5 


Nd II 


4825, 


.48 


0. 


.18 


-0.880 


5 


Nd II 


4832, 


.27 


0. 


.56 


-0.920 


5 


Sm II 


4777, 


.85 


0. 


.04 


-1.340 


5 


Sm II 


4815, 


.80 


0. 


.18 


-0.770 


5 



'^(l) McWilliam (1998); (2) Bard, Barisciano, & 
Cowley (1996);(3) Thevenin (1989);(4) Meggers, 
Corliss, & Scribner (1975); (5) VALD data base 
(Piskunov et al. 1995) ;(s) solar gf derived. 

^^Hyperfine structure included according to 
McWilliam (1998). 
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Table 6. PREDICTED C/0 RATIO AND / VALUE IN EXTRINSIC AGB STARS 





1.5 Mq 




[Fe/H] 


(C/0)™ 


/ 


+0.00 


1.45 


1.60 


-0.12 


1.65 


1.15 


-0.30 


2.00 


0.79 


-0.40 


4.28 


0.25 


-0.52 


5.02 


0.21 


-0.60 


5.55 


0.18 


-0.70 


6.28 


0.16 


-0.82 


7.36 


0.14 


-1.00 


9.18 


0.11 


-1.30 


13.29 


0.07 


-1.60 


33.08 


0.03 
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Fig. 1. — Best synthetic fit (solid line) to the observed (filled circles) spectrum of the star 
Z Psc in the range A4800 - 4830 A. The main atomic absorption lines are indicated. Note 
that there are still many spectral features not well reproduced by our synthetic spectrum. 

Fig. 2. — As Figure 1 for the star in the region of the Ba II line at A4934 A for different Ba 
abundances: no Ba (upper dashed line), 12 + log (Ba/H)=2.13 and 3.1 (best fit, sohd line). 

Fig. 3. — Distribution of R = F{12fim)/ F {2.2 fim) for Tc-ycs (open histogram) and Tc-no 
(shaded histogram) N stars. The Tc-no stars with R > 0.6 are collected in the last bin. 

Fig. 4. — The (K-[12], K-[25]) color-color diagram, where open circles stand for Tc-yes and 
filled circles for Tc-no N stars. Filled diamonds are stars in our sample with no information 
about Tc (whether from this work or the literature). Filled squares are the carbon stars of 
J-type studied by Abia & Isern (2000). J-stars do not show Tc. The solid line represents 
black body colors, with crosses corresponding to temperatures 4000, 3000, 2000 and 1500 
K. The error bars represent a typical uncertainty of 5 and 10% on the 12 and 25 fxin fluxes, 
respectively. Note that there is not any segregation between different groups of carbon stars 
in this diagram. 

Fig. 5. — Comparison of the observed mean low-mass s-element (Y,Zr) [Is/M] enhancement 
against metallicity with theoretical predictions for a 1.5 Mq TP-AGB star. The theoretical 
predictions shown are for 0-enhanced stellar models at C/0— 1. The upper and lower curves 
(dashed fines) fimit the region allowed by the models according to the different ^^C-pocket 
choices (see text). Note that several stars coincide in the same data point. The difficult stars 
are plotted with smaller symbols (see text). 

Fig. 6. — Same as Figure 5 for the observed mean high-mass s-element (Ba,La,Nd) [hs/M] 
enhancement as compared with theoretical predictions for a 1.5 Mq TP-AGB star. 

Fig. 7. — Same as Figure 5 for the observed mean [hs/ls] ratio (signature of the neutron 
exposure) as compared with theoretical predictions for a 1.5 Mq TP-AGB star. 

Fig. 8. — Observed trend of the Is element abundances [Is/M] vs. [hs/ls] for Galactic disk 
intrinsic and extrinsinc AGB stars. The solid curves refer to envelope models with different 
s-process efficiency, as monitored by the A^(^'^C)/N(^^Fc) ratio (here normalized to the case 
that fits the main componet in the solar system). The dashed lines connect the points corre- 
sponding to the fourth and eighth dredge-up episode to make clear how the stars distribute 
along the TP-AGB evolutionary sequence. Again, our difficult N stars have been plotted 
with smaller symbols. 

Fig. 9. — A detailed reproduction of the observed abundances in the carbon stars S Set, AQ 
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And, EL Aur (left panel), UV Aql, Z Psc, and UU Aur (right panel). See text for details. 
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H 1 h 
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El Aur 
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40 60 

Nuclear Charge Z 



80 
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0.5 



2 

1.5 h 

Cl. 

\, 1 h 

0.5 P 
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El. 



Ed 



1.5 

1 h 

0.5 7 

- 



-0.5 



T 1 1 1 1 

UV Aql 

model: [Pe/H] = 0. M^l.SM^, ST 




Z Psc 

model: [Fte/H] = 0, M=1.5Ma. ST* 1.3 



4 
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H 1 h 





H 1 ^ 



UU Aur 

model: [Fe/H] = 0. M=1.5Mq. ST/ 1.5 
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I 
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40 60 

Nuclear Charge Z 
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Table 4. DEPENDENCE OF DERIVED ABUNDANCES ON MODEL ATMOSPHERE PARAMETERS 





ATeff = ±200 K 


Alog= ±0.5 


= ±0.5 kms-i 


AC/0 = ±0.05 


Ai2C/i^C= ±10 


ACNO/H = ±0.2 


A[M/H]= ±0.3 


[Ba/H] 


±0.02 


±0.15 


TO.05 


±0.10 


±0.0 


±0.01 


±0.05 


[La/H] 


±0.0 


±0.20 


TO.30 


±0.10 


±0.0 


±0.01 


±0.10 


[Ce/H] 


±0.05 


±0.30 


TO.40 


±0.10 


±0.0 


TO.10 


±0.10 


[Nd/H] 


±0.02 


±0.20 


TO.40 


±0.10 


±0.0 


TO.10 


±0.05 


[Sm/H] 


±0.0 


±0.20 


TO.40 


±0.0 


±0.0 


±0.01 


±0.10 
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Table 5. ABUNDANCES DERIVED IN PROGRAM STARS 



Star [M/H] [Rb/M] [Sr/M] [Y/M] [Zr/M] [Ba/M] [La/M] [Ce/M] [Nd/M] [Sm/M] Tc 



AQ And 


0.02 


0.3 


0.6 


0.7 


0.5 


0.1 


0.1 


0.3 


0.0 no 


AWCyg 


0.0 


0.2 


0.4 


0.3 


0.5 


0.0 


0.3 


<0.5 


yes 


EL Aur 


-0.06 


0.4 


0.5 


0.7 


0.8 


0.2 


0.2 


0.1 0.2 


no 


HK Lyr 


-0.10 


0.2 


0.3 


0.8 


0.7 


0.8 


0.8 


0.7 


yes 


IRC -10397 


0.10 


0.1 


0.3 












no 


lY Hya 


-0.80 


0.5 


0.6 


<0.7 


<0.7 


<0.5 








LQ Cyg 


0.25 


0.2 


<0.5 


0.8 


0.7 


0.4 


0.4 


0.5 


<0.4 no 


RVCyg 


0.0 




0.0 


0.3 




0.0^ 






no 



2 



Table 5 — Continued 



Star [M/H] [Rb/M] [Sr/M] [Y/M] [Zr/M] [Ba/M] [La/M] [Ce/M] [Nd/M] [Sm/M] Tc 



RX Set 


-0.05 


0.4 


0.4 


0.7 


0.8 


0.7 


0.7 




0.8 


<0.5 


no 


S Set 


0.01 


0.5 


0.8 


0.7 


0.5 


0.2 


0.1 




0.2 




yes 


SS Vir 


0.0 




0.0 


0.5 


0.4 


0.3 


0.4 




0.3 




doubf 


SY Per 


-0.30 


0.1 


0.5 
















doubf 


SZ Sgr 


-0.04 


0.1 


0.4 


0.9 


0.8 


0.8 


0.9 




0.9 


<0.6 


doubf 


TT Cyg 


-0.1 






1.2 


1.0 


0.8 


0.9 




0.9 


<0.6 


doubf 


TY Oph 


0.10 


0.2 


0.6 


0.7 


0.6 


0.6 


0.5 




0.6 




doubf 


U Cam 


-0.09 






0.8 


0.5 


0.5^ 


0.5 


0.2 


0.5 


0.3 


no 



3 



Table 5 — Continued 



Star [M/H] [Rb/M] [Sr/M] [Y/M] [Zr/M] [Ba/M] [La/M] [Ce/M] [Nd/M] [Sm/M] Tc 



U Hya 


-0.05 


0.5 


0.8 


1.3 


1.1 


1.1 


0.9 


0.6 


<1.0 


<0.7 


yes 


UU Aur 


0.06 


0.4 


0.4 


0.6 


0.6 


0.3 


0.2 




0.3 


<0.2 


no 


UV Aql 


-0.06 


0.6 


0.5 


0.7 


0.9 


0.6 


0.6 




0.7 




yes 


UX Dra 


-0.20 


0.0 


0.8 


0.8 


0.7 


0.7 


0.4 




0.4 




yes 


V Aql 


-0.05 


0.3 


0.5 






0.2^ 










no 


V CrB 


-1.35 


-0.2 


0.3 


0.4 














no 


V Oph 


0.0 


-0.2 


0.2 


0.3 


0.4 


0.1 


0.1 


-0.05 


0.2 






VX Gem 


-0.15 




0.6 


0.5 


0.6 


0.8 




<1.0 




<1.0 


yes 
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Table 5 — Continued 



Star [M/H] [Rb/M] [Sr/M] [Y/M] [Zr/M] [Ba/M] [La/M] [Ce/M] [Nd/M] [Sm/M] Tc 



VY Uma 


-0.1 


0.3 


0.5 


0.8 


1.1 




1.0 




0.8 


<0.6 


yes 


V460 Cyg 


-0.04 


0.4 


0.5 


0.7 


0.8 


0.8 


0.7 




0.8 


<0.4 


yes 


V781 Sgr 


0.10 


0.1 


0.5 


0.3 


0.3 


0.1 


0.2 




0.1 




yes 


W CMa 


0.0 


0.2 


0.6 


















W Ori 


0.05 




0.0 


0.3 


0.2 


0.3 


0.3 




0.1 




no 


Z Psc 


-0.01 


0.6 


0.9 


1.0 


1.0 


1.0 


1.1 


0.6 


0.9 


0.8 


yes 



^The stars RV Cyg, U Cam and V Aql show weak Merill-Sanford (SiC2) band absorptions not included in the analysis 
that could affect the derivation of the Ba abundance. 



